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ABSTRACT 



A reduced form of the patch conic method has been employed to 
determine the initial orbital elements of a particle launched or 
ejected from the Moon's surface with any arbitrary starting condi- 
tions. The reduction was obtained by considering the selenocentric 
velocity asymptotes. Explicit and tractable analytic functions have 
been derived for the geocentric and Jacobi energies, angular momentum, 
standard orbital elements, and conditions for Mbon-to-Earth trajec- 
tories. Percents of randomly ejected material which initially strike 
earth, are in retrograde orbits, or go into heliocentric orbits have 
been obtained. The results are compared with results obtained by a 
numerical integration program for several different situations . 
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INTRODUCTION 

The problem of the time history of motion in space of trajectories 
originating on the Moon's surface has received attention in recent 
years. There are three fundamental areas in which the problem arises 
which are currently under investigation. The Moon has been considered 
by several authors as the origin of tektites (Ref. l) ; a lunar source 
was suggested by Whipple as a possible origin of the dust cloud about 
the Earth (Ref. 2) j and the space program has given rise to the prob- 
lem of returning manned or unmanned spacecraft from the Moon's surface. 

The validity of the idea that there must be some kind of large 
chunks of material from the Moon on Earth, and some amounts of lunar 
dust arriving daily on Earth has been established by studies of crater- 
forming processes. Several large craters on Earth have been shown to 
be generated by meteoritic impact. (Ref. 3 ). Shoemaker's analysis 
(Ref. 4) of the ejecta pattern of the lunar crater Copernicus sup- 
ports the hypothesis that it is of impact origin and shows that some 
of the energy of impact is consumed in the formation of secondary 
impact craters caused by ejecta from the primary crater which have 
acquired large kinetic energies. It is plausible to assume that some 
fragments acquired kinetic energy in excess of the Moon's gravita- 
tional potential energy. Further, Gault et al's experimental studies 
(Ref. 5 ) indicate that bombardments by smaller meteors may cause a 
small fraction of the ejected material to acquire escape velocity from 
the Moon. Therefore it is of interest to determine the likelihood of 
this material arriving on Earth, and the general behavior of all 
material ejected into space. 
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The usual method of investigating Moon-to-Earth trajectories is 
to employ the patch-conic method to establish preliminary orbits, and 
then, if precise results are required, to refine the orbits by 
utilizing numerical integration methods. The patch-conic method, which 
was apparently originated by Tisserand and first applied to problems 
involving space travel to the Moon by Egorov (Ref. 6), considers the 
particle to be in a Keplerian orbit about the body exerting the major 
force on the particle. 

The sphere of influence is the surface about a planet or moon 

within which that body exerts a stronger force than the perturbative 

force of any other body on the particle. The surface is nearly a 

m 2/5 

sphere with radius given by r = fi (-) where m/M < 1 and ft is the 
distance between m and M. For m = mass of the Moon and M = mass of the 
Earth, r ~ 66,000 km. Similarly, a sphere of radius 1 million km 
around the earth describes the region in which the earth's gravita- 
tional field predominates over that of the Sun. 

In the case of Moon-to-Earth trajectories, the patch-conic method 

considers that the particle is in selenocentric Keplerian hyperbola 

» 

from the Moon's surface to the edge of the sphere of influence. Then 
the position and velocity coordinates with respect to the Moon are 
converted into geocentric position and velocity coordinates, taking 
into account the position and velocity coordinates of the Moon 
relative to the earth. The Earth's gravitational field is then con- 
sidered to be the sole deteiminant of the consequent motion. 

The purpose of this paper is to establish the nature of the initial 
geocentric orbits of particles leaving the Moon's surface with arbitrary 
initial conditions. In order to accomplish this, a limiting form of 
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the patch-conic method is developed into analytic forms describing 
the geocentric motion of any particle leaving the Moon’s surface. 

The formulas establish the approximate initial orbital elements; more 
accurate work or consideration of motion over a long duration requires 
numerical integration procedures involving all of the relevant 
perturbations. However, the approximate formulas provide physical 
insight, preliminary Moon-to-Earth trajectories, the average percentage 
of material ejected from the Moon that arrives on Earth within a few 
days, the distribution of retrograde versus direct geocentric orbits 
and the percentage of material with geocentric kinetic energies in 
excess of the energy stored in the gravitatio nal fields of the Earth- 
Moon system. 



FORMULAS 


A. CONCEPTS 

The object of this paper is to perform a survey of the important 
characteristics of the motion of bodies leaving the surface of the 
Moon after an initial, instantaneous thrust. The discussion will 
center particularly around the problem of dust ejected by meteoritic 
bombardment of the Moon's surface. There are a wide variety of initial 
conditions which may be imparted to particles ejected from the Moon 
during crater formation. Craters are observed to cover the surface 
of the Moon in an irregular manner. Each individual crater sprays 
out material at a wide range of angles and velocities. 

There are six coordinates necessary to describe the initial con- 
ditions of the particle at ejection from the Moon. The position coor- 
dinates may be given by the radial distance from the Moon's center, 
lunar latitude and longitude; the velocity coordinates are given by 
the speed, the azimuth direction and the elevation angle. Of these 
coordinates, only the radial distance from the Moon's center may 
reasonably be taken to be a constant, equal to the radius of the Moon. 
There remain five variables to be studied for an indefinite number of 
craters. By using some reasonable assumptions, it will be shown 
that the major features of resulting geocentric conditions can be 
obtained with only three explicit variables. 

The first simplification is to consider that material is ejected 
from the Moon uniformly on a surface of a fictitious sphere surrounding 
the Moon. The trajectories arising from any particular crater will 
not be considered. Instead, it will be assumed that, on the average, 
the location of points where the ejecta pierce the sphere is randomly 
distributed. 
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If the fictitious sphere is at a large enough distance from the 
Moon so that the motion is near the asymptotic state of hyperbolic 
trajectories, the velocity vectors will not be randomly oriented in 
space. The vectors will be clustered around a trajectory where the 
velocity and position are aligned in the radial direction away from 
the Moon, that is, the trajectories on the average appear as spokes 
centered around the Moon (Fig.l). This assumption permits the velo- 
city vector to be defined with the same angular coordinates as the 
position vector, reducing the number of variables from five to three. 

The three variables which play an important part in the analysis 
are v r , 0, ep. The magnitude of the residual velocity v r is defined 
to be 


2 

v 0 


2 

V 

esc 


( 1 ) 


2 2GM m 

where v Q is the ejection velocity v = is the classical 

m 

two body velocity of escape from the Moon's surface. The residual 

— 4 

velocity v^ makes an angle 0 with the velocity vector of the Moon 
v m< The angle 0 is taken to be 0° when v^ is in the direction of the 
Moon's motion and l80° when opposite to the Moon's motion. In spheri- 
cal coordinate terminology 0 is the "polar angle", although the Moon's 
physical pole is perpendicular to the pole of this velocity 
coordinate system. The angle cp specifies the azimuthal component of 
v^ measured in a cone formed by 0 = constant. The plane of the Moon's 
orbit around the earth is used for 0 to establish the value cp = 0. In 
this approximation, no distinction will be made between the earthward 
side and back side, or between the northern and southern hemispheres; 
therefore 0° £ cp < 90°. (There is an important difference between the 
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forward and backward velocity hemispheres which is accounted for by 

taking 0 to be between 0° and 180° . ) 

The components of the position and velocity vectors are to be 

obtained with respect to an "inertial" Earth-centered coordinate 

system as a function of v , 0, cp. First the components are expressed 

relative to a Moon-centered system rotating around the Earth (Fig. 2). 

Choosing £ to be in the direction of the Moon's angular momentum 

vector with respect to the Earth, i to be pointing to the Earth, 5 

to be in the direction of the Moon's antapex, that is, opposite to 

the Moon's velocity vector, and remembering that the pole of the 0, cp 

coordinate system is in the -3 direction, the residual velocity vector 
— ♦ 

v f may be written 

v r = f v r sin0coscp - 3 cos0 + £ v r sin0sincp (2) 


The position vector with respect to the Moon is assumed to be vector 
of zero amplitude because of the paradoxical- sounding assumption that 
the particle is an infinite distance from the Moon but at lunar dis- 
tance from the Earth. In order to obtain coordinates with respect 
to the Earth-centered inertial system, the axes are translated to 
the Earth, maintaining the orientation of the axes so that there is 
now an instantaneous geocentric position vector 


= -{R 


m 


( 3 ) 


In the inertial coordinate system, the Moon is moving at an instan- 
taneous rate described by 


v = 
m 


3 v. 


m 


00 
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The vector addition of the particle's velocity with respect to the 
Moon and the Moon's velocity with respect to the Earth 




provides the geocentric particle velocity v as 

& 


a function of v , 
r* 


e, cp 


— * A A 

v = i v sinScoscp - j (v cos© + v ) 
g r T ° r m 

+ £ v^ sinQsincp 

The foregoing conditions may be viewed as describing a fictitious 

physical situation in which the Moon has "vanished", removing from the 

particle all of the potential which is due to the lunar graviational 

field. This is a reduced form of the patch conic technique. The 

more general form of the patch conic analysis could be established by 

o o 

restoring terms involving R^, ©, cp. 

The definition of v r (Eq. 1) provides a lower limit on v Q in order 
to have particles orbiting in Earth -Moon space. 

v ° a v esc * 2 -37 

This minimum value of v Q will be used herein; however, it may be 
pointed out that the use of v defined from two body considerations 

Q SC 

is somewhat misleading. Particles with ejection velocities much less 

than v will describe ballistic trajectories which strike the Moon 
esc 

after a brief flight. As v 0 increases, the particles will describe 
incomplete Keplerian ellipses. As these ellipses become larger and 



more elongated, they will be increasingly subject to Earth perturbations. 
In some cases, perilune may be raised enough to permit the particle 
to become a temporary lunar satellite. When v Q becomes large enough 
to satisfy the condition 


2 

Vo * 


esc 



the particle has sufficient energy to cross the Moon’s sphere df 

influence of radius R^, and may do so, depending upon the value of 

the angular momentum. Between, this value of v 0 , and the minimum 

value of Vo obtained above, there is a "fuzzy" region -where it can’t 

be definitely said whether the particles "escape" the Moon's sphere 

of influence or not. Furthermore, above v 0 = 2.37, 'the particles 

cannot be said to "escape" the Moon in the two body sense that they 

can never return to the Moon. The term "escape velocity" will be 

2 GM 

used to mean that v Q = -p— . 


The major effect of the Moon on the initial trajectories has been 
accounted for by removing the gravitational potential and the Moon’s 
geocentric velocity from the initial conditions of the particle. The 
remaining geocentric coordinates of the particle will be operated on 
with two-body (particle and earth) formulas to establish the basic 
types of initial trajectories. The most important single parameter 
of an orbit is the energy, so that calculation will be presented first. 

The law of conservation of energy provides the two -body formula 


, 2 
E = hr 


ri £ 2 


r 2 


(9) 


where E is a constant and the subscripts refer to sets of coordinates, 
Obtaining from (6) 

222 2 

Vi = v - v + 2 v v cos 9 + v 
g r r m m 

and ri = R , the value of E can be obtained 
nr 


( 10 ) 


2 2 GM e 

E = iv + V V cos 0 + 4v - - 

c r r m c m R 

m 


( 11 ) 


In the two body model, the particle is said to have escaped when 

E > 0. If E < 0, the particle will trace a Keplerian ellipse or circle. 

GM 

2 e 

Substituting v = = — , the condition E > 0 can be stated as 
m R 7 


m 


2 2 

v + 2 v v cos0 - v ^0 
r r m m 


( 12 ) 


9 



If cos 8 = +1, this condition requires that v ^ (/2 - l) v for the 
particle to escape, or, if cos = -1, then v r ^ 1/2 + l) for the 

particle to escape. Therefore if v^ < .4 v ffi , no particles will 
escape, if .4v m < < 2.4 v , the particle will escape if the 

direction of v^, is sufficiently close to the Moon's direction of mo- 
tion; if v r > 2.4 v , all particles will escape. The term "escape" 
is used here to mean that the particle has sufficient kinetic energy 
to overcome the Earth's gravitational potential and recede to "infinity" 
during the initial orbit . 

The variation of geocentric energy for a given value of ejection 
velocity, depending on the value of 0, may seem to be a violation of 
the conservation of energy. Indeed, it is shown below that the absolute 
or Jacobi energy depends only on the ejection velocity. However, the 
immediate geocentric energy depends upon the vector addition of the 
residual velocity vector with the Moon's velocity vector. This is 
illustrated schematically in Fig. 3- 

A more realistic escape condition may be derived by considering 
the system of Earth, particle and Sun. If the particle crosses the 
earth's sphere of influence (~10 6 km), the particle will be in es- 
sentially heliocentric orbit with a theoretical possibility of re- 
turning to the Earth-Moon system at some later time. The energy re- 
quired to leave the Earth-Moon system is somewhat less than to es- 
cape to infinity (Eq. ll) 


e = yi 


GM 


+ V V COS0 

r m 


K - ^ 


m 


ft 

h 


The values of eacape to infinity given in Eq. 11 or from the 
Earth-Moon system given in Eq. 13 do not provide information as 


to the possibility of ultimately escaping the Earth-Moon system. The 
use of the Jacobi integral of energy obtained from the restricted 
three-body theory applied to a system assumed to consist of the Earth, 
Moon and particle provides more insight. Expressed in the non-rotating 
coordinates system, and assuming the barycenter to be located at the 
Earth’s center of mass, the Jacobi integral is 


„ GM GM 

I v 2 - n' Cxy - yi) = ^ J - C 

e m 


v 

where n’ = — , r a and r m are the distance of the particle to the 


m 


m 


Earth and Moon, respectively, and C is the constant of integration. 

Under the previous assumptions discussed earlier in this paper, the 

GM 

particle is far enough from the Moon so that = 0 while remaining 

m 

at nearly the same distance from the Moon as the Earth so that 

2 2 

r = R . Then v = v is given in Eq. (10) and x = -R , so the 
e m g m* 

expression for the Jacobi constant reduces to 


i 2 i 2 _ 

* V r " * v m “ R 


GM 


m 


A consistent calculation is obtained by considering the Jacobi 
integral in the rotating system of coordinates: 

GM 


A - 4 n ,2 R 2 = 
c m 


■e v.. 


R 

m 


- C 


which obtains again Eq. 15. 

Comparing the Jacobi constant with the expression for the two- 


body energy yields 


Since E is not expected to be a true constant of motion, •while C is, 

the changes' in E must be contained in the term v v cos9. While it 

r m 

is not within the scope of the present paper to discuss the long-term 

behavior of the particles, a situation may be imagined in order to 

clarify the meaning of Eq. 17 . Consider a particle which initially 

has E < 0 but v > {/, 2 - l)v . If at some later time the particle 
r m 

passes through the Moon's sphere of influence and the Earth's gravita- 
tion has a negligible effect on the particle during this passage, 
then, by the laws of two-body mechanics, the residual velocity magni- 
tude after passage must equal the residual velocity magnitude before 
passage. However, the particle's trajectory will bend around the Moon, 
the amount of bending depending on the closeness of the encounter with 
the Moon. Therefore the value of cos 9 will change, the value of E 
may fluctuate within a range of 2 v^v^, hence a particle which 
originally has a value of E < 0 may acquire enough additional energy 
by a passage through the Moon's sphere of influence to bring E above 
zero. That is, a particle tracing perturbed Keplerian ellipses may 
become hyperbolic with respect to the Earth-Moon system by this 
mechanism. So E may be said to define the limits of a particle 's 
motion during the early phases of its history, while C provides 
information on the ultimate limits. 
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c. 


ANGULAR MOMENTUM 


The geocentric angular momentum is a useful quantity to obtain, 
especially in order to evaluate the probability of a particle coming 
directly to earth. The angular momentum (for a particle of unit mass) 
is 

A = Rxv = j R v sin8 sincp + k R (v cos0 + v ) 
g J m r T mr m 

The condition necessary to establish a trajectory which intersects 

the center of the earth is that there is zero angular momentum. This 

condition is met when the 3 and Is components above are set equal to 

zero: 

cp = 0° | 

-V COS0 = V J 

r m 

The orbit of a particle with those initial conditions will be a straight 

line traveling toward or away from the earth lying in the plane of the 

Moon's orbit. For a retrograde orbit, the ^-component must be negative, 

that is v cos0 + v <0. 
r m 

The square of the magnitude may be obtained from Eq. 18 

A 2 = |R x v g j 2 = Rj [v^(l-sin 2 0 cos 2 cp) + 2v r v m cos0 + /] 

This function, although dependent on the three variables v , 0, cp, is 

still tractable enough to be used for analytic investigations. For any 

value of cp, the condition for minimum angular momentum may be found by 
2 

differentiating A with respect to 0, 

2 2 

= -2v cos 0 - 2 v v sin0 
r r m 

2 2 
+2v sin0 cos0 sin cp 
r 


1_ d(A 2 ) 

R 2 d0 
m 


( 18 ) 


(19) 


( 20 ) 


( 21 ) 
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Setting this derivative equal to zero provides the value of 0 where 
2 

A (or A) is a minimum for fixed values of cp: 

-v 

v cos0 = =— 

r 2 

cos cp 

If cp = 0, the formula ( 20 ) is obtained. The two equations cp = 0 

and v cos 0 = -v provide the minimum of the function 
r m 

A(v r , 0, cp) = 0. 

The magnitude of the angular momentum vector can be used to pro- 
vide the initial conditions of trajectories which intersect the Earth. 
The Tnnvi m um value of ang ula r momentum Ag a particle can have and be 
able to strike the Earth is given by 


ISe x " *E t s s1q 9°° " ^ 

where R is the radius of the Earth, and v is the velocity of the 
e s 

particle at the Earth's surface. The magnitude of v g can be computed 
from the conservation of energy: 




Setting 2GM (5- - J-) = vf - 122.9 km 2 /sec 2 and using v (v ,0,cp) 

"e a 6 


we obtain 


Comparing A^ with A (Eq. 20), the values of 9 which establish an orbit 
which just grazes the Earth's surface are expressed as a function of 
v and cp: 


COS0 = 


-v ( 1 -a 
m 

2 

V COS cp 


J: 


v 2 (l- 
m 


a? 2 ) + cos 2 cp [o^v 2 - (l-cv 2 ) (v 2 


v 2 ) 

nr 


2 2 
V COS 

r 


v cos 
r 


«P 


«p] 

(27) 


2 — 

where o = = 2.75 x 10 4 is small so the first term is roughly 

m 

equal to the condition (22) for minimum angular momentum.' Calling the 
first term of cos0 f Q and the second term f, it may be said that those 
trajectories that have values of 0 such that 


f° + f £ cos0 ;> f° - f (28) 

will strike the earth during the initial orbit. Incidentally, cos0 is 
necessarily negative in this region. This region is enclosed between 
the dotted lines in Tig. 4 for the case when cp = 0. As cp increases 
from zero, the dotted lines move closer together until finally the 
region vanishes altogether. The maximum value of cp which will permit 
and Earth-strike is 


coscp 


max 


-b + s/b 2 -4-ac 


2a 




where 


2 

a = v r 

2 2 , 2 2 2 . 
b = O' V & - (V + V ) (1-0' ) 

2 / i 2.2 

C = V l l-cv ) 
m 

Equation (29) is obtained by setting the radical in (27) 


0. 


(29) 

(29a) 

(29b) 

(29c) 
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There is a discontinuity in Eq. (27) at cp = TT/2., of the order 
2 

of (TT/2) . An alternate form of cos0 may be used when cp = TT/2 


cos0 = - 


2 

V - 

m 


2 2 
a v 

a 


2v v 
r m 


taking 


a 

l-o 


= a 


This equation is only needed when, as computed from (30), |cos9) < 1 
or, solving for the roots of the equality. 


( 30 ) 


v - a v < v ■< v + a v . (31) 

m arm a 

This is the region of v r on Fig. 4 where the first dotted line has 
appeared and the second has not. 
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D . ORBITAL ELEMENTS 


The geocentric orbital elements a , e,. i, iu, Q, v are provided 
below as a function of the coordinates v , 9 , cp. Although the semi- 
major axis a and the eccentricity are more complicated functions 
than those for the equivalent constants of "energy" and "angular 
momentum" described in the previous two sections, they are provided 
for completeness. The inclination i establishes the conditions 
which cause retrograde geocentric orbits . The right ascension of the 
ascending node is constant within the accuracy of this analysis, a 
surprising result which is verified by numerical integration. The 
argument of perigee u) is the negative of the true anomaly v as a 
consequence of the assumptions in this analysis. 

Using equations (3) and (6) which provide r and v in tenns of 

§ 

v r > cp, the usual formulas converting cartesian position and 
velocity vectors into orbital elements (Ref. 7) may be applied. The 
reduction GM = v^ R^ is used throughout this section. The semimajor 
axis, given by 


a = 


GM 

2 GM - rv" 


becomes 


a = - R 


m 


m 2 0 2 

v + 2v v oos0 - v 
r r m m 


The eccentricity, which is 


e = [ 1 - 


]Rxv 
GM a 


2 * 


( 32 ) 


( 33 ) 


( 34 ) 
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3 


<M • 


2 2 t* 2 

( 1-sin 0 cos cp) + 2 — — - cos© (2 - sin 0 eos 2 cp) 


V r /, £ 2 2 

+ — — l4cos 0 + sin 0 cos c 


The angular elements i and fi are obtained from the direction components 
of the unit angular momentum vector 


R X v 


is x \\ 


From Reference 7, the expressions 


li = sin fJ sin i 
x 

-fi = cos n sin i 

y 

fi = cos i 
z 

may be used to obtain the form 


tan i = ± 


(h 2 + h 2 ) 
x y ' 


which provides a simple expression for the inclination as a function 
of v r , 0, cp 


tan i = 


v sin 0 sin cp 
r 

V cos 9 + v 
r m 


If only the positive value of the square root is taken, the denomina- 
tor provides the sign of the fraction. If v cos e + v <0, the 

r m 
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orbit is retrograde as indicated previously in Eq. 18. 

The right ascension of the ascending node Q given by 


h 

tan Q = - — (4o) 

y 

is zero for all values of v , 0, <p since the 1 component of angular 
momentum is zero. This feature is a consequence of the coordinate 
system. Physically speaking the node lies along the Earth-Moon line, 
or the x-axis of the system. If some other reference point in the 
Earth-Moon system is used to define the zero position of the node, then 
the node will be a non-zero constant. 

The true anomaly is obtained from 

tan v . I g y - glf 8 ■ — = f 0 *) 

R x v] -GMr 


where 

s = v sin0 coscp fv (l-sin 2 0 eos 2 cp) + 2v v cos0- +v } 
r “ r T r m m 

and 

22 2 2 2 

t = v cos 0 + 2v v cos0 + v sin 0 sin cp 
r r m r T 


(42) 

(43) 


Unfortunately, these complicated functions are apparently irreducible. 
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The argument of perigee is related in general to the true anomaly 
by the formula 

id = u-v (44) 

where 

sin u = (4-5) 

z 

Since z =0, the argument of perigee is equal to the negative of the 
true anomaly. 

This analysis has provided four new independent variables 
a, e, i, a) and two dependent or constant quantities v and Q. 

Normally, six independent variables are required to specify particle 
motion. That there are only four is a consequence of utilizing only 
four independent variables v ,• 0, cp and R^, or, expressed physically, 
assuming that the radial and velocity vectors with respect to the 
lunar system are coincident. (Although R^ has the mathematical 
property of being an independent variable, in most of this discussion 
it is only necessary to consider one value) . 
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RESULTS 


A. INTRODUCTION 

The Earth-Moon system provides an interesting application of the 
preceding theory because the velocity of escape v g from Hie Moon's 
surface is larger than the Moon's geocentric orbit velocity v . 

Except for the Earth's Moon, the escape velocity from the natural 
satellites of the solar system is less than the orbital speed of the 
satellite about the planet. (See Table I). In this case, the 
residual velocity of an ejected particle is likely to be small compared 
with the orbital velocity of the body from -which escape occurred, so 
that the particle's orbit about the planet will be similar to that of 
the satellite . 

For a particle to leave the Moon, the ejection velocity must be 
much larger than the orbital speed. Small changes in the ejection 
velocity imply relatively larger changes in the residual velocity 
vector. Therefore the residual speed of a particle ejected from the 
Moon may be a sizeable fraction of, or equal to, or greater than the 
Moon's geocentric orbital speed. Performing the vector addition 
(Eq. 5) of a residual velocity which is comparable in magnitude but 
different in direction from the Moon's velocity yields a geocentric 
velocity which will be quite unlike the Moon's velocity vector. If a 
set of arbitrary residual velocity vectors is considered, it will be 
found that the net geocentric velocities will vary widely among them- 
selves. When transformed into geocentric orbits, a variety of types 
will be obtained, often varying sharply with small changes in the 
speed v r or the direction (expressed with 0 and tp) of the residual 
velocity vector. 
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The geocentric orbits obtained will be illustrated in this section 
in several ways. The resulting geocentric orbits may be categorized or 
described in several ways. The following pages discuss the results in 
three frameworks suitable for different applications. The energy of the 
orbits which will exist for one or more orbits in the earth-moon system. 
The location of and characteristics of regions where the initial 
trajectories intersect the Earth are useful for many problems, including 
manned lunar return missions, and the origin of tektites. The general 
characteristics of the orbital elements obtained previously may be of 
interest. 


TABLE I 

Ratio of Escape Velocity to Orbital Velocities of the Massive Natural 
Satellites of the Solar System 


Planet 

Satellite 

Escape Velocity 

Orbital Velocity 

Ratio 

Earth 

Moon 

2.37 km/sec 

1.02 km/sec 

2.3 

Jupiter 

Io 

2.28 

17.31 

.13 


Europa 

1.99 

13.73 

.15 


Ganymede 

2.83 

IO.87 

.26 


Callisto 

2.23 

8.20 

.27 

Saturn 

Titan 

2.99 

5.57 

.54 

Neptune 

• Triton 

2.80 

4.4o 

.64 
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B. ENERGY 

The energy of an orbit is perhaps the most important single 
describer of an orbit. If the energy is small, the particle is bound 
to the Earth-Moon system and can never escape. If the energy is. high 
enough, the particle will leave the Earth-Moon system during its first 
orbit and will not fall back, although it may later pass through the 
Earth -Moon system during a chance encounter. 

Previously, two different "energies" were obtained (Eq. 9 and 
Eq. 12) . The Jacobi integral of energy C is the true indication (with- 
in the limits of the restricted three -body theory) of the ma-xi muin pos- 
sible distance from the Earth -Moon system after a long period of time. 
The Jacobi constant depends in this analysis, only upon v , that is, 

only upon the energy acquired during the ejection process. The value 

GM 

of the Jacobi (expressed in canonical units where ^ = l) as a fune- 

m 

tion of v^ is given in Fig. 5. It should be noted that the definition 
of the Jacobi constant C is inverted with respect to the definition 
of energy. That is, the larger the velocity becomes, the smaller C 
becomes, whereas the two-body energy becomes larger with increasing 
velocity. 

Michael (Ref. 8 ) shows the curves of zero velocity around the L 4 
and L 5 for the Earth-Moon system (his definition of C contains a fac- 
tor of two compared with C as used in this paper) . The value of C at 

the L 4 is approximately the value of C arrived at when v =0: 

r 

however, because of the mode of analysis involved in obtaining C in 
this paper, it should not be concluded that a particle with zero 
residual velocity will become trapped in the L 4 position. Indeed, this 
analysis makes no distinction in potential along the Moon's radius. 
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It is instead assumed that the Moon has "vanished ", removing all of its 
gravitational potential from the particle’s energy. The utility of the 
Jacobi constant obtained herein is, recognizing that the particle 
may reencounter the Moon after one or more orbits, in indicating how 
the two body energy may be modified at some later time. 

The immediate, or two-body (Earth and particle) energy E depends 
upon two parameters, v and 0. Because of this, values of E may be 
plotted (Fig. 4) as a function of v r and 0. Since the probability that 
a spherically symmetric distribution will provide a given value of 0 
is proportional to cos0, the 0-axis is marked in cosine units. The v r 
axis is linear with respect to the ejection velocity v 0 . A better 
choice would be a velocity axis proportional to the distribution of 
ejection velocities as obtained from cratering theory or experiments, 
but that is not within the scope of this paper. However, it may be 
recognized that the lower ejection velocities are more probable. Since 
the patch conic technique has been found to give poorer results at low 
velocities (Ref. 9), the plot begins at 2.k km/sec. 

The escape condition (Eq. 12, E = 0) is shovn as Line 1 on Fig. 4. 
At the lower ejection velocities, most of the particles exiting from 
the rear hemisphere (with respect to the Moon’s motion) do not have 
enough energy to "escape to infinity". As the ejection velocity in- 
creases, the area above Line 1 steadily shrinks until v Q = J.k km/sec. 
it vanishes altogether. 

A less stringent condition than "escape to infinity" may also be 
of interest. Then energy required to leave the Earth-Moon system 
(Eq. 13) is obtained for values of v r and 0 along Line 2. Having suf- 
ficient energy to leave the region where the Earth’s gravitational field 
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predominates and go into essentially heliocentric orbit does not neces- 
sarily imply the particle will. The angular momentum must be low as 
well. However, Line 2 is useful in defining an upper limit to the 
region where particles can immediately escape in any sense. 

The percentages of orbits which are elliptic with respect to the 
Earth may be obtained by integrating 
TT/s 180 

i dcp J sin0 d9 where 

0 0 . 

01 

2 2 

v - v 

01 = cos" 1 is obtained from Eq. 12 and 

r m 

normalizing by dividing by the area of a spherical quadrant. This 
1+COS01 

obtains — — X 100% which is plotted Fig. 6 as a function of v Q . 

A similar procedure yields the percentage which has less energy than 
required to leave the Earth -Moon system during the first orbit. 

The angular momentum is particularly useful in demonstrating the 
conditions necessary for a particle to come directly to the Earth. 
Furthermore it indicates whether an orbit is direct or retrograde. 

It was shown (Eq. 39) that the condition -v^cos© = v ffl represents 
the demarkation line between retrograde and direct orbits. This con- 
dition is met along Line 3 in Fig. 4. Above line 3 the orbits are 
retrograde, below they are direct. Four classes of orbits may be 
distinguished: 

I) direct and hyperbolic, the predominate group 

II) direct and temporarily trapped 

III) retrograde and temporarily trapped 

IV) retrograde and hyperbolic. 

The percentages of retrograde orbits as a function of ejection velocity 
is given in Fig. 6. 


It was established also that the equation -v cos 9 = v is also 

r m 

one of the two conditions establishing a trajectory which intersects 
the center of the Earth. If it is assumed that Fig. 4 is specialized 
to represent only the cp = 0° plane, then Line 3 provides the initial 
conditions which yield zero angular momentum with respect to the earth. 
The region where particles may pass below the Earth's radius on the 
first orbit is between the dotted lines 4 and 5 shown on Fig. 4, pro- 
vided that cp=0. The orbits between lines 3 and 5 will be retrograde; 
between lines 3 and 4 direct. 
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C . MOON-TO-EARTH TRAJECTORIES 

The problem of Moon-to-Earth trajectories has been studied 
previously by authors interested in manned space flight and in the 
origin of tektites. In the former category, Penzo (Ref. 9) and 
Dallas (Ref. 10 ), working mainly with the patch conic technique, iden- 
tified many of the major characteristics of such trajectories. In the 
latter category, Varsavsky (Ref. 11), in a paper several years ago, 
found several Moon-to-Earth trajectories by numerical integration. 
However, his results were too meager to provide much insight into the 
problem. It is of interest to apply the formulation obtained earlier 
in this paper to obtain general statements about such trajectories. 

Trajectories arriving on Earth during the first orbit were shown 
in the preceding section when cp = 0°, that is, the particle moves 
entirely in the plane of the Moon's orbit around the Earth. The con- 
ditions necessary to obtain a Moon-to-Earth trajectory can be illus- 
trated better by utilizing a coordinate plane containing 0 and cp. The 
axes of 0 and cp are formed by projecting the rear velocity hemisphere 
onto. a polar coordinate plane about a pole coincident with the Moon's 
antapex. The radial spokes correspond to values of cp; the pole is at 
9 = l80°, the outermost circle is 0 = 90°. The cp = 90° line divides 
the hemisphere into l/k sphere facing the Earth, and l/4 sphere facing 
away from the Earth . 

By using Eq. 27, the sets of v , 0, cp which provide direct access 
from the Moon to the Earth can be obtained. By considering a particular 
value of the ejection velocity v Q , curves enclosing small regions in 
the 0 , cp plane are established within which the 0, cp points will yield 
Moon-to-Earth trajectories for that value of v D . These regions, which 
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will be called 'Velocity Strike Zones" in this paper, are plotted in 
Fig. 7 for several specific values of v . 

Several major features may be pointed out : 

1. A minimum value of ejection velocity somewhat higher than 
escape velocity from the Moon is required to establish 
these trajectories. It is 



( 47 ) 


or, approximately, 

■ v ' s= v - a v (48) 

r m a 

which requires v^s .8^4 km/sec, or v Q ^ 2.51 km/sec. 

2. When v Q is low, but above the minimum, there is one velocity 
strike zone surrounding 0 = l80°. As v 0 becomes larger, 
the zone begins to contract in the middle and elongate along 
the 0-axis, becoming somewhat dumbbell shaped. Next, the 
zone fissions into two separate zones, one on either side of, 
but not including, 0 = l80°. The zones are oval shaped, and, 
in this approximation, are identical to each other. The 
trajectories stemming from the zones on the earthward side 
of the hemisphere are headed toward perigee; those on the 
outer side pass through apogee before coming to perigee. 

As v 0 becomes larger, the ovals move further outward along 
the «p = 0° axis, shrink slightly in size, and become more 
circular. When v 0 is large enough to cause a hyperbolic 
trajectory with respect to the Earth- Mson system, the zones 
on the hemisphere away from the Earth vanish, as the 
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particles will never pass through the perigee phase of their 
orbit. As v 0 becomes very large (compared with v ) the 
zone approaches 6 =90° and shrinks to the limit imposed by 
geometrical size of the Earth. This limit can be obtained 
by considering Eq. 2J when the conditions v » v 
(v ~ 11 km/sec) hold. The equation becomes cos9 =• 0 ± or, 
that is, the center of the zone lies along the Earth -Moon 
line and the zone extends on either side an angular distance 
n = Rg/R^. This is the geometrical limit which was used by 
Urey (Ref .12) to obtain a negligibly small percent for the 
amount of ejected material arriving on earth. 

The velocity strike zones lie close to the plane of the 
moon's orbit. This is in agreement with Penzo r s observa- 
tion that the asymptotic velocity vectors lie within 10° 
of the Earth-Moon plane. 

The boundary of the velocity strike zone was defined by 
considering the maximum value of angular momentum necessary 
to strike the Earth or the Earth r s atmosphere. This condi- 
tion provides trajectories which just graze the Earth's 
atmosphere, i.e., have shallow entry angles. Within this 
boundary lie concentric curves which yield constant values 
of angular momentum. The values of angular momentum become 
lower as the curve becomes closer to the center of the zone. 
At the center, the angular momentum goes to zero as dis- 
cussed earlier. Therefore, a desired value of entry angle 

can be maintained by moving along a curve generated by 

Rg cosp 

replacing a = with a 1 = — ' • • where cos$ is the 


angle of entry into the atmosphere with respect to the 
horizon. 

Most of the results presented above have been obtained previously 
by authors using the usual form of the patch conic method. Dallas 
gives a, discussion of the properties of the loci of exit points from 
the sphere of influence for a given reentry angle of Earth, which are 
analogous to the curves obtained from Eq. 27. Penzo has also pointed 
out many of the above features. The advantage of the simplified version 
of the patch conic method described in this paper is that it provides 
elementary analytic equations to work with. 

The choice of the v , 0, cp coordinate system has some disadvantages 
for this part of the problem. The discontinuity at cp = TT/2, while man- 
ageable, causes clumsy results — for instance, the maximum value of cp 
for Earth strikes drops suddenly from TT/2 to about 1^°. However, the 
use of the conventional latitude and longitude would result in an 
algebraically more complicated formula . 

The percent of ejecta which arrives on Earth during the first 
orbit can readily be obtained as a function of ejection velocity by 
ecomparing the area of the Velocity Strike Zone with the area of the 
unit sphere. Because of symmetry, it is only necessary to consider 
the portion of the unit sphere bounded 0° & 0 £ 180, 0° < cp < 90°. 

The area of this quarter-sphere is Aj/ = TT. The area of the i 

/ 4 ^ 

portion of the velocity strike zone for a given value v r contained 
in this quarter sphere is given by 


cp 0 2 

T max 


_ p XUfcLA p 

I Vsz = J J slne de dc P 


0 0 ! 

where 0i = cos" 1 (f° + f) and 0 2 = cos -1 (f° - f) are defined in Eqs. 


C^9) 


50 



( 50 ) 


(27) and (28) and cp is given in Eq. 29. 

max 

The total percent is given by 

P = X 100# 

a 7 * 

The computational procedure varies, depending on the value of v^: 

1. If f° + i' H < -1, or equivalently, v < v - or v 

i Hi El 

(v 0 = 2.51 lcm/sec) , then no particles will arrive on 
Earth and P 0. 

2. If f° + f > -1 but f° - f < 1, then 

TT/2 - e 

P = j (f° + f + 1) dtp (51) 

where £ is a small number which is used to avoid the 

discontinuity in Eq. (27) at TT/2. Since the discontinuity 
2 

is of order 7^-, the numerical accuracy can be maintained 
close enough to TT/2 so that little contribution to P is 
lost. 

3 . When f° - f“ > 1, v > v + of v a (v Q = 2.66 km/sec) then 
algebraic cancellation can occur and 


P = 


100 

~ir 


c 


max 


f dcp 


Although there is an abrupt change in the upper limit at 

v = v + a v . when cp drops sharply from TT/2 to about 
r m b. max 

15°, P is unaffected and varies smoothly thru the transi- 
tion point. 


(52) 
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The results of the computations described above are shown in 
Fig. 8 . From the minimum velocity of 2.51 km/sec required for an 
earth strike, the percentage rises sharply to the maximum value. 

An ejection velocity of ~ 2.6 km/sec will send the maximum amount of 
material directly to the earth, about 3 - 3 $. The curve decays some- 
what slower, to provide a range in velocity of about .3 km/sec which 
will send at least 1$ of the ejecta directly to the Earth. At about 
v Q = 2.96 km/sec, there is a discontinuous drop in the percent because 
the geocentric orbits are now hyperbolic and only those trajectories 
passing through the Velocity Strike Zones facing the Earth will ever 
go through perigee. After this the curve approaches very gradually 
the asymptotic limit imposed by the geometric size of the Earth. 

These results are in good agreement with Chapman's (Ref. 13) , 
which were obtained by computing over 4000 trajectories with numerical 
integration methods. He presents in graphical form the percentage of 
ejecta arriving directly on Earth as a function of ejection velocity. 
The percentages are given for material ejected at positions on the 
moon where ten large craters are located. The curves have the same 
shape and location of the peak as in Fig. 8 . Because individual 
craters have greater or lesser amounts of material which will pass thru 
the Velocity StrikeZones, the peak-values of the percent curves vary 
between 1 and 6 $. It can be seen, however, that the average peak would 
be about 3 . 5 $, as is indicated in this paper. 
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D. ORBITAL ELEMENTS 

Results have been obtained for the geocentric orbital elements 
using the analytic formulas developed in Eq. 32-45. Typical examples 
of the variation of the orbital elements are presented graphically 
in this section. 

The variations of the semimajor axis with 0 are shown in Fig. 9 
for specified values of v . Only the positive values of a are shown 
in this figure. Eq. 33 shows that a(0) is a conic, and the curves 
look hyperbolic. There is a flat bottom approaching 0 = l80°, with 
a sharp rise going to the left at the lower values of 0. At the point 
indicated by Line 1 on Fig. 4 the semimajor axis becomes infinite and 
at still lower values of 0 it recedes from minus infinity with nega- 
tive values. As the ejection velocity increases, the curves become 
steeper. 

The eccentricity as a function of 0 is shown in Fig. 10 for 
specified values of the ejection velocity at cp = 0° and in Fig. 11 
for specified values of cp at v = 2.6 km/sec. Values of e < 1 are 
elliptic orbits, e > 1 are hyperbolic and, e = 1 parabolic. It can 
be seen from the formula for perigee height, q = a(l-e), that the 
condition that e be near one is a necessary but not a sufficient 
condition for an Earth-strike. In the case of v Q = 2.6 km/sec the val- 
ues of 0 near l80° where e ~ 1 yield Moon- to -Earth trajectories, while 
the value of 0 near 90° where e ~ 1 do not. 

The inclination is shown in Fig. 12 as a function of cp for 
several values of v 0 and 0. The sharpest rises occur near the Velocity 
StrikeZones and, as will be discussed in the next section, may in 
fact be even sharper than indicated by the analytic formula for inclina- 
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tion. It is interesting to note that the high inclinations relative 
to the plane of the Moon's orbit, both retrograde hnd direct, are 
possible and are especially apt to occur in the case of Moon-to-Earth 
trajectories . 

The argument of perigee is shown as a function of 0 in Fig. 1J for 
v 0 = 2.6 km/sec. The curve varies between 0° and 180° because of the 
assumed symmetry of the front side and back side of the Moon. However, 
negative values should be given to the argument of perigee if the 
asymptote passes through the rear side of the Moon's sphere of influence, 
and positive when through the Earth side. 
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E. ACCURACY 


The patch conic method has been successfully utilized in the study 
of trajectories for lunar missions, including manned space flight and 
needs no justification here. Empirical corrections have been utilized 
to improve the accuracy of the patch conic result (Ref. 9 ); but because 
the aim of this paper was to provide the maximum possible simplicity in 
computation, and because specific answers such as landing point or time 
* of flight were not required, no attempt was made for high accuracy. 
Indeed, complicated refinements would destroy the utility of the 
present work. If improved accuracy is needed, the usual form of the 
patch conic, or if necessary numerical integration, should be re- 
sorted to. 

Still it is of value to provide some justification as to the 
general correctness of the formulas obtained within this paper. 

For this purpose, use was made of a numerical integration program 
(ITEM Ref. 14) of good precision. ITEM is an interplanetary trajec- 
tory flight program which is based on the Encke method of integrating 
trajectories influenced by several gravitational or other perturbing 
forces. It has been applied to investigate or predict the motions of 
several artificial satellites as well as lunar missions. 

Several situations were compared between the two programs. 

Usually spoke orbits, as were assumed in developing the analytic 
theory, were computed until they had left the Moon’s sphere of influence 
(after about 10 hours flight time). The geocentric orbital elements 
provided at this point were compared with the elements obtained from 
the analytic theory. 
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A series of runs was made along the plane of the Moon's equator 
to indicate the changes in the orbital elements with changes in 0 . 

The values of v = 2.5 km/sec, 0 every 30° and cp = 0° were used. Fig. 14 
shows the comparative results for perigee. The ZCEM values were taken 
after 48 hours of flight. (A curve constructed after 18 hours of flight 
time deviated a little more from the curve constructed by the analytic 
method, but not enough to illustrate clearly on a graph of this scale.) 
Two curves are obtained from ITEM to test the effect of the symmetry 
of the orbits leaving from the backside with those leaving from the 
Earth side of the Moon. The curves are nearly congruent and have the 
same shape so the symmetry hypothesis is reasonable for this order of 
accuracy. The semimajor axis is shown in Fig. 15 (negative values are 
for hyperbolic cases). The agreement is satisfactory. 

The velocity of 2.5 km/sec, being just above the velocity of 
escape from the Moon, represents the region where the patch conic method 
is poorest. (Ref. 9). The same set of runs using an initial velocity 
of 3.05 km/sec showed deviations between the three curves -which were 
too small to represent on graph paper of the same scale as Fig. 14. 

The region near Earth-strike was investigated with values of 0 oc- 
curing every 2° between 118° and l4o°. (Fig. 16). The two programs 
yield similar curves with a slight lateral displacement which may be 
due to the differences in the Moon's distance. (Variations in the 
locations of the Velocity Strike Zones have been observed to occur with 
changes in the Moon's position. (Ref. 9 ). 

A general test was made by considering spoke orbits with v Q = 2.8 
km/sec at 0 = 90° and <p = 0-90 every 15 ? (When 0 = 90 , the longitude 
is 180° and the latitude is equal to cp) . The comparative results are 
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NOTE: First row is figures obtained from ITEM; second row from analytic method. 


indicated in Table II. Here the semimajor axis is negative, indicating 
a hyperbolic orbit. There is some divergence between the two programs 
in the semimajor axis, but ITEM finds it to be a constant to within 1$. 
The eccentricity, inclination and argument of perigee are correct to 
about three places, but the true anomaly divergences increasingly as 
the initial conditions are moved upwards from the plane of the equa- 
tor from the negative value of the argument of perigee. 

The relative elements of Moon-to-Earth trajectories as obtained 
by ITEM and by the analytic formulas were compared. Trajectories were 
computed for 0 = 140°, 9=0°- 13° and vE = 2.7, 2.8 km/sec. This 
region is within the Velocity Strike Zones for v = 2.7, and 2.8 km/sec. 
Therefore the results of this set of runs gives some indication of the 
accuracy of Moon-to-Earth trajectories. Table III shows the values of 
the six orbital elements, as obtained from numerical integrations (ITEM) 
of the spoke orbit and by the formulas discussed earlier in this paper. 
The orbital elements which vary along a perturbed trajectory, were 
the osculating elements ten hours after launch when the particle was 
about 60,000 km from the Moon. A launch time on April 6, 1965 at 
2?: 00 UT was arbitrarily assumed in order to initiate the ITEM 
computation. The semimajor axis differs between the two programs by 
about 1000 km; however, the theoretical prediction that the semimajor 
axis is independent of, changes in 9 provided that a constant value of 
0 is maintained is verified within 9 The eccentricity agrees to 
four places. The inclination disagrees badly, apparently because the 
denominator in Eq. 39 goes to zero in the case of a Moon-to-Earth 
trajectory. Then terms involving the z-coordinate, which were omitted 
in the derivation of Eq. 39, may become important. The difficulties 

are compounded by the necessity of taking the arctangent of a large 
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number . The node is remarkably constant as predicted; the change of 
node with change of. velocity is because the ITEM program happens to 
use the instantaneous Earth-Moon line as a reference. The true 
anomaly is within 1° of the negative of the argument of perigee and 
these angles are within 1° when comparing between the two programs. 
Because the true anomaly varies continuously during the trajectory 
and because a situation cannot be obtained in a numerical integra- 
tion program which is comparable to the assumption that "the moon has 
vanished instantaneously" absolute agreement with the prediction 
that uu = — v cannot be expected. 

For completeness, a study was made of a set of orbits which were 
not spoke orbits but more realistic orbits, physically speaking. 

The spoke conditions 0 = 135°, cp = 45° and v Q = 2.7 Im/sec were chosen 
and, assuming an inclination of 45° to the reference plane, starting 
conditions on the Moon’s surface were calculated for several assumed 
angles of ejection to the Moon’s horizon. The trajectories were 
numerically integrated past the edge of the Moon's sphere of influence. 
The lunar coordinates of the particle did not correspond exactly to 
the assumed spoke conditions. These deviations may he ascribed to 
the influence of Earth perturbations during selenocentric phase of the 
trajectory and to the fact that the asymptotic conditions with respect 
to the Moon have not been reached exactly, or to loss of place accuracy 
in the hand computation Obtaining the surface conditions, or to the 
difference between the Moon’s orbital plane and its equatorial plane. 

Table IV presents the results of the computations. Since the 
latitude at exit from the sphere of influence, which would be 30 ° 
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TABLE IV 

COMPARISON OF VARIOUS NON-SPOKE ORBITS AT ONE POINT IN 0, cc SPACE 

Ejection Velocity = 2.7 km/sec 
Spoke Conditions: 6 = 135° cp = 45° 


•ITEM Computations: 

Selenocentric Conditions at Lunar Surface: Selenocentric Conditions 10 hours (53000km) a 

Ejection 


Elevation Angle 

Latitude 

Longitude 

Velocity 

Elevation Angle 

Latitude 

0.0 

-44930 

-1719700 

1.35 

86?09 

27953 

30.0 

-24.54 

- 120.880 

1.35 

86.65 

29936 

45.0 

- 11.12 

- 101.880 

1-35 

87.25 

28.95 

60.0 

+ 2.92 

- 87 .O 80 

1.35 

88.02 

28.88 

Analytic Method: 






90.0 



1.28 

90.00 


Ejection Velocity 

= 2.8 km/sec 



at 9 hours 


0.0 

-44.00 

-164.87 

1.54 

86.5 

28.6 

30.0 

- 22.23 

-120.05 

1.34 

87 .O 

29.1 

45.0 

- 10.12 

-100.46 

1.54 

87.5 

28.7 

60.0 

+ 3.75 

- 86.26 

1.54 

88.2 

28.7 

Analytic Method: 






90.0 



1.48 




Geocentric Orbjjtal Elements 


a 

e 

3P2092 

. 742914 

375924 

.711023 

390152 

•735950 

384515 

.741741 

322373 

.715246 

520954 

.744389 

485696 

.718152 

525564 

.757075 

519872 

.759296 


i 

G 

84.896 

-6.127256 

85 . 544 

-6. 18^060 

84.126 

-6 . 183294 
-6 . 154400 

1 

84.755 

80.30 

1 

95.105 

-6.036072 

96.857 

-6.076914 

95.071 

-6.067245 

95.535 

-6.029534 


U> 

V 

-138.642 

142.044 

-136.917 

140.589 

-136.942 

140.585 

-138.192 

141 . 842 

-145.63 

145.63 

-122.473 

126.398 

-123.049 

127.115 

-123.713 

127.707 

-124.576 

128.596 


4-09842 


.708167 


92.37 


-131.42 


131.42 



for 9 = 135 ° , cp = ^5°, was as much as 2 l/2° off, some error is 
introduced in comparing the exit elements with the analytic elements. 
Under these conditions, the deviations between the methods are in the 
second place. This is sufficient accuracy for some purposes, and with 
care, more accuracy could be obtained. 

There are several sources of deviations between HEM and the 
analytic method. The geocentric orbital elements provided by ITEM 
are not at all constant but vary continuously along the trajectory. 

First, the energy of the particle at the edge of the Moon's sphere of 
influence still contains a perceptible component of the Moon's poten- 
tial energy. The accuracy might be increased by matching elements at 
a later time in ITEM, or by removing only the difference in potential 
energy with respect to the Moon. Second, the eccentricity of the Moon's 
orbit induces changes in R m and v . Third, the perturbations by the 
Earth during the Moon phase of the trajectory are neglected. Fourth, 
the particle is given an x-coordinate only so that terms depending on 
the y or z coordinate are neglected. Fifth, the particle is assumed 
to be moving along a hyperbolic asymptote with respect to the Moon, an 
assumption that improves as the ejection velocity increases. 

In summary, the analytic method provides a good approximation to 
the geocentric orbital elements over a wide range of initial conditions. 
In order to evaluate precisely how much accuracy is provided by the 
method, precise specifications must be made as to what is being evaluated 
and when . The accuracies indicated by the previous tables and graphs 
are sufficient for the purposes of this paper. The accuracy of the 
method can probably be improved for particular purposes. 


CONCLUSION 


A good description of all possible trajectories arising from 
launch or ejection from the Moon's surface has been obtained in the 
form of a few elementary algebraic formulas . It has been shown that 
velocity asymptotes in the instantaneous direction of the Moon's 
motion will yield orbits that are hyperbolic with respect to the 
Earth-Moon system even if the ejection velocity is barely above that 
necessary for escape from the moon. It was indicated that, as the 
ejection velocity increases, there is an enlarging area centered about 
the Moon's apex where the velocity asymptotes will provide hyperbolic 
geocentric orbits. When the ejection velocity is large enough, all 
geocentric orbits will be hyperbolic. 

The demarkation between retrograde and direct geocentric orbits 
was given. The condition establishing trajectories with zero angular 
momentum was obtained. 

The initial conditions for Moon-to-Earth trajectories has been 
established analytically. The geometric conditions for a given 
initial velocity have been computed and shown graphically. The formulas 
obtained are in good agreement with numerical results described by 
previous authors. The minimum velocity (v 0 = 2.55 km/ sec) for a Moon- 
to-Earth trajectory requires that the velocity asymptotes lie near the 
Moon's antapex. The Velocity Strike Zone encompassing the largest area 
is about v Q = 2.65 km/sec. While these trajectories are too long in 
duration to be practical for manned space flight, they might be 
desirable for unmanned spacecraft returning from the Moon. 

The significance of "gravitational focusing" on Moon-to-Earth 
trajectories has been discussed in the literature. The foregoing 
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analysis shows that the relative velocity of the Moon about the Earth 
is of more importance in determining the characteristics of Moon-to- 
Earth trajectories. .Although the amount of the velocity added by 
falling through the Earth's gravitational field is large (~ 11 km/sec), 
its effect oh Eq. 27 is small, since it is multiplied by a and the 
product is small (~ .06 km/sec). 

The initial geocentric orbital elements were found to be in rea- 
sonable agreement with results of numerical integrations. Improvement 
could probably be made by allowing for the eccentricity of the Moon’s 
orbit, or with other corrections indicated herein. 

The general conclusions made on the basis of analysis were sub- 
stantiated to a remarkable degree by the numerical integration runs. 

1) The semimajor axis is independent of cp. This was true to 
four places in Table II. 

2) The right ascension of the ascending node is zero (or con- 
stant, depending on the reference point) . This prediction is confirmed 
to almost as many places as are carried in the numerical integration 
program. (The change of Cl with change in velocity is a spurious ef- 
fect due to the particular coordinate system utilized in ITEM. ) 

3) The prediction that the initial true anomaly is equal to the 
negative of the argument of perigee is well verified. This is a 
geometrical effect occuring because the particle starts in the 
reference plane. Since the particle has traveled for l/2 a day be- 
fore the Earth-referenced true anomaly is computed, there is a dif- 
ference of about a degree in v and -u). 

This form of analysis is useful because of its functional sim- 
plicity and probably could be extended for applications to more closely 



defined, problems . The method can probably be used to generate initial 
conditions for the patch-conic iterations to find Moon-to-Earth 
trajectories and to provide analytical partial derivatives to home in 
on the final trajectory. With a modification to allow for the effect 
of arbitrary perilune, the method might be applied to the problem of 
circumlunar trajectories. 
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